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Abstract 



The process e + e~ — > e + e~Z/7* is studied with the OPAL detector at LEP at a 
centre of mass energy of y/s = 189 GeV. The cross-section times the branching ratio 
of the Z/7* decaying into hadrons is measured within Lorentz invariant kinematic 
limits to be (1.2 ±0.3 ±0.1) pb for invariant masses of the hadronic system between 
5 GeV and 60 GeV and (0.7 ±0.2 ±0.1) pb for hadronic masses above 60 GeV. The 
differential cross-sections of the Mandelstam variables s, t, and u are measured and 
compared with the predictions from the Monte Carlo generators grc4f and PYTHIA. 
From this, based on a factorisation ansatz, the total and differential cross-sections 
for the subprocess e7 — ► eZ/7* are derived. 
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1 Introduction 



In this paper the reaction e + e~ — > e + e~Z/7* is studied using the OPAL detector at LEP 
and the cross-section times branching ratio for the decay of Z/7* into hadrons, denoted 
as cr ce , is measured. In this reaction a quasi-real photon is radiated from one of the beam 
electrons and scatters off the other electron producing a Z/7* as shown in Figure [l| This 
process was measured for the first time |1[ with the OPAL detector. The observable final 
state, (e)eff , consists of the scattered electron, e, and a fermion pair, ff , from the Z/7* 
decay while the other electron, (e), usually remains unobserved in the beam pipe due to 
the small momentum transfer squared, \p\ 2 , of the quasi-real photon. 

From the cross-section a ee the cross-section of the subprocess e7 — > eZ/7*Q denoted by 
<T e7 , is determined. This is the first measurement of the cross-section a crf (\fs) for values of 
v§ equal to or greater than the Z-mass. The process e7 — > eZ/7* is the same as ordinary 
Compton scattering with the outgoing real photon replaced by a virtual photon 7* or a 
Z. 

The cross-section a ee is given by the convolution of the cross-section <t C7 with the 
photon flux D n (z, s) 

a ee (s) = / dz D n (z,s) da n (s) 
•j [1 



1 • 




where z = s/s. 

For Z boson or 7* production in Compton scattering e{k)~f(p) — > e(k') Z/7* (p') of real 
photons (p 2 = 0), the cross-section depends on the Mandelstam variables s = (k + p) 2 = 
(p> + k') 2 , t = (k' - k) 2 = (p' - p) 2 and u = (p' - k) 2 = (k' - p) 2 [§ 

2 4d + i). (2) 

us uj 

The variable m q q is the invariant mass of the quark-anti-quark pair the Z/7* decays into 
and for m q q = the well known terms for ordinary Compton scattering remain. 

A singularity at u = is introduced by the virtual electron propagator in Figure |l](b) 
as the typical transverse momentum scale of the scattered Z/7* bosons is small . For 
incoming quasi-real photons (p 2 ~ 0) in ep or e + e~ collisions, the dominant regulating 
effect for this divergence is not the electron mass, but small, non-zero, incoming photon 
masses squared p 2 . Via the replacement [0, |J 

2 

- 1 2 m qq 2 fn\ 

u — » u + p -^- — m e (6) 

in the denominator of Equation Q, both the photon mass and the electron mass are 
included in the propagator. 



1 Charge conjugation is implied throughout the paper except when otherwise stated. 



A simple equivalent photon approximation (EPA) 



LUO,s)epa 



'cy\ 



a 1 + (1 - z) 
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m e 



(4) 



where the integration is performed over the small photon virtualities, leads to an effective 
on-shell incoming photon flux. This overestimates the cross-section by a factor of two || . 
The p 2 spectrum of the incoming photons either has to be retained fully or modified to 
describe the process properly. The modified EPA, denoted by EPA is given by ||: 



D ei [z, u) E p A 



a 1 + (1 
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(5) 



In any case, the results will be sensitive to the modelling of the p 2 spectrum. In this 
paper the theoretical expectations are represented by Monte Carlo event generators using 
different approaches for obtaining the p 2 spectrum of the incoming photons. These are 
compared with the experimental data. The comparisons include the distributions of the 
Mandelstam variables s, t, and u as well as other characteristic variables, like m q q, and 
E e , the energy of the scattered electron. 

After giving a description of the data used for this analysis and of the OPAL detector, 
a signal definition is given. Using kinematic invariants a part of the cross-section a ee 
is defined as signal. Thereafter the selection of the signal events is described and the 
total cross section a ee within the signal definition is calculated. Using the same selection 
differential cross-sections dcr ee and d<j e7 are determined using an unfolding technique. 



2 Data and detector description 



The analysis uses 174.7±0.2 (stat.)±0.3 (syst.) pb _1 of data collected during 1998 with the 
OPAL detector at LEP at a centre of mass energy of y/s ~ 189 GeV. A detailed description 
of the OPAL detector may be found elsewhere || and only a short description is given 
here. The central detector consists of a system of tracking chambers providing charged 
particle tracking over 96% of the full solid angleQ inside a 0.435 T uniform magnetic field 
parallel to the beam axis. It is composed of a two-layer silicon microstrip vertex detector, 
a high precision drift chamber, a large volume jet chamber and a set of z chambers 
measuring the track coordinates along the beam direction. A lead-glass electromagnetic 
(EM) calorimeter located outside the magnet coil covers the full azimuthal range with 
excellent hermeticity in the polar angle range of | cos#| < 0.82 for the barrel region and 
0.81 < | cos 6*| < 0.984 for the endcap region. The magnet return yoke is instrumented 
for hadron calorimetry and consists of barrel and endcap sections along with pole tip 
detectors that together cover the region | cos#| < 0.99. Four layers of muon chambers 

2 The OPAL coordinate system is defined so that the z axis is in the direction of the electron beam, 
the x axis points towards the centre of the LEP ring, and 9 and <f> are the polar and azimuthal angles, 
defined relative to the +z- and +cc-axes, respectively. The radial coordinate is denoted as r. 



cover the outside of the hadron calorimeter. Electromagnetic calorimeters close to the 
beam axis complete the geometrical acceptance down to 24 mrad, except for the regions 
where a tungsten shield is designed to protect the detectors from synchrotron radiation. 
These calorimeters include the forward detectors (FD) which are lead-scintillator sandwich 
calorimeters and, at smaller angles, silicon tungsten calorimeters |7| located on both sides 
of the interaction point. The gap between the endcap EM calorimeter and the FD is 
instrumented with an additional lead-scintillator electromagnetic calorimeter, called the 
gamma-catcher. The tile endcap || scintillator arrays are located at 0.81 < | cos 6\ < 0.955 
behind the pressure bell and in front of the endcap ECAL. Four layers of scintillating 
tiles || are installed at each end of the detector and cover the range of 0.976 < | cos#| < 
0.999. 



3 Signal definition and event simulation 



3.1 Signal definition 

The predominant signature of the signal process in the final state (e)eqq is one electron, 
two hadronic jets from the Z/7* decay and large missing momentum in the direction of 
the beam pipe due to the escaping electron. The cross-section is peaked at low \u\ where 
the scattering angle of the signal electron is large, i.e. in the backward^ direction, and 
its energy is small. Furthermore, \u\ — > implies that the Z/7* is emitted close to the 
forward direction. As a consequence a huge part of the cross-section lies outside of the 
acceptance of the OPAL detector; therefore the signal is defined within kinematic limits. 

The process e + e~ — > e + e~Z/7* and its subprocess e7 — > eZ/7* can also be measured 
at a future e + e~ linear collider || ||. There this process will be the dominant source of 
real Z production. Furthermore this subprocess can be observed in ep collisions []2], |I~0| 



where the beam proton emits a bremsstrahlung photon. The relevant quantity for the 
e7 collision is Vs, the centre of mass energy in the e7 rest-frame. In order to be able to 
compare the results of this analysis with the measurements of other experiments Lorentz 
invariant quantities are used for the definition of the signal. This is in contrast to [|I|] 
where the signal was defined by the geometrical acceptance of the detector rather than 
by Lorentz invariant quantities. Consequently the results from the previous paper cannot 
be compared directly with the ones presented here. 

The signal is defined by the two diagrams shown in Figure [l] within additional kine- 
matic limits as detailed below. Further processes like the ones shown in Figure ^] leading 
to an (e)eqq final state are treated as background even if they satisfy our signal definition. 

The Feynman diagrams for t-channel Bhabha scattering with initial or final state 
radiation are identical to the e + e~ — > 7*e + e~ Compton scattering diagrams. While the 



3 The forward direction is denned by the initial direction of the electron radiating the Z/7* 



Bhabha events with initial-state radiation of a virtual photon correspond to the w-channel 
7*ee diagram, Bhabha events with final state radiation are equivalent to the s-channel 
7*ee events. The cross-section for Bhabha scattering diverges for p 2 — > 0. This divergence 
is regulated by a finite p' 2 of the radiated 7* but it still causes the cross-section to be 
largely peaked at small \t\ = \(p' — p) 2 \. Bhabha scattering with 7* radiation may best 
be characterised by two energetic electrons (small p 2 of the exchanged photon) and a 
preferably low- momentum 7* (small p' 2 ) in the e + e~ centre of mass system, leading to 
small \t\. In the observable phase space of the e + e~ — > e + e~Z/7* process on the other hand, 
the energies of the incoming photon (p 2 ) and outgoing Z/7* (p' 2 ) are sizeable and their 
momenta prefer opposite directions, leading to large negative values of t. We therefore 
require the absolute value of the kinematic invariant \t\ to be larger than 500 GeV 2 to 
define our signal. 

The square of the four-momentum transfer of the quasi-real photon, \p 2 \, is required 
to be less than 10 GeV 2 to ensure that the electron emitting the quasi-real photon stays 
within the beam-pipe. As \p 2 \ is usually small, this requirement does not reduce the cross- 
section by much. In order for the EPA from Equation |5| to provide correct results the 
virtuality of the quasi-real photon needs to be the smallest virtuality in the process. This 
is guaranteed by requiring \u\ of the electron in Figure [jj)) to be larger than 10 GeV 2 , the 
cut value on \p 2 \. This cut mainly rejects events which would be very difficult to select 
because either the energy of the scattered electron is small or the scattering angle is very 
close to the beam direction. 

In order to avoid the region of hadronic resonances with all its uncertainties in the 
simulation of the spectrum we require the square of the invariant mass of the Z/7*, m 2 



qq; 



to be greater than 25 GeV 2 . The kinematic limits for the signal are summarised in table |l|. 



Angle and energy of the signal electron: 


\i\ > 500 GeV 2 


Mass square of the quasi-real photon: 


\p 2 \ < 10 GeV" 


Virtuality of electron: 


\u\ > 10 GeV" 


Mass square of the qq system: 


m 2 q > 25 GeV" 



Table 1: Cuts used for the definition of the signal 
The interdependence of the Mandelstam variables is given by 

t = -I( S -m 2 q )(l-cosr) (6) 

u = -I(,- m 2 _)(l + cosr) (7) 

where cos#* is the scattering angle of the Z/7* with respect to the e7 axis in the e7 
rest-frame. Defining the kinematic region of the signal within |t| > 500 GeV 2 and «iqq> 
25 GeV 2 is an effective cut on the centre of mass energy in the e7 rest-frame at Vs > 
22.9 GeV. The kinematic invariant t is completely determined by the four-momentum of 
the electron and u is determined by the hadronic decay products of the Z/7*. Neglecting 



the mass of the electron one obtains 

i = -2 EE C (l + Q c cos# c 



u -2EE m (l -I^QeCOsflqd+mJi, (9) 



E 

-^qq 



with E being the energy of the beam electrons, E e , cos9 e and Q e the energy, scattering 
angle and the sign of the charge of the electron, E qq and p qq the energy and momentum 
of the hadronic system. The cut on t is therefore a cut in the [E e , — Q e cos# e ] plane as 
depicted in Figure 0. Since cos# c > —1 there is a hard cutoff on E c > 1.3 GeV. 



3.2 Signal Simulation 

For the generation of the e + e~ — > e + e~Z/7* signal events two different Monte Carlo gen- 
erators, grc4f [IT] and PYTHIA 6.133 [I2~|, are used. 



The grc4f Monte Carlo generator is linked to GRACE, an automatic Feynman dia- 
gram computation program. The total and differential cross-sections are obtained from a 
phase space integration of the matrix element, which is calculated from all diagrams cor- 
responding to a given initial and final state. All fermion masses are non zero and helicity 
information is propagated down to the final state particles. Also a subset of diagrams 
can be chosen. Here only diagrams according to the signal definition have been used. A 
sample of events corresponding to about 30 times the data luminosity has been analysed. 

In PYTHIA the cross-section is calculated according to Equation |2] including the 

regularisation given in Equation |3| to avoid the divergency in the matrix element. In 

contrast to grc4f the matrix element for the process e7 — > eZ/7* and the modified EPA as 

given in Equation |5| are being used. In PYTHIA a cutoff on p t , the transverse momentum 

of the Z/7* with respect to the axis of motion of the electron and the photon in the e7 

rest-frame, is applied. The default cutoff of 1 GeV has been removed in order to include 

the full phase space. This has been made possible by introducing the new regularisation of 

Equation [3] into PYTHIA 6.133. A further replacement is made to ensure the cross-section 

does not become negative: 

2 

i^i-^-i (10) 

™qq 

A sample of events corresponding to approximately 11 times the data luminosity has been 
used. 

For both Monte Carlo generators parton showers and hadronization of the final quarks 



are performed by JETSET |T2|] with parameters tuned to the OPAL data [|13| . 



Within the kinematic limits defined above the cross-section a ee is predicted by grc4f 
to be (1.77 ± 0.02) pb, while the corresponding value from PYTHIA is (1.92 ± 0.03) pb. 
The errors are statistical only. 

8 



Figure |] shows the distribution of p t , wi q q, E e and t on generator level for the two 
Monte Carlo samples after applying the signal definition. PYTHIA predicts a higher cross- 
section, mainly at small electron energies and large values oi pt. In the m qq - distribution 
the contributions from the 7* and the Z are well separated. 



3.3 Background Simulation 

The main contribution to the background comes from two-photon hadronic processes, 
e + e~ — ► e + e~+ hadrons. These events are divided into three subsets, depending on the 
virtualitiesQ, q 2 and p 2 , of the photons and consequently the number of beam electrons 
observed ("tagged") in the detector. For the low momentum transfer processes ("un- 
tagged"), both q 2 and p 2 are small; these are simulated using PYTHIA 5.7. Where the 
momentum transfer of one of the photons is large ("single tagged"), i.e. q 2 is large and p 2 



small, HERWIG [|14| is used. The PHOJET |L5| generator is used for the processes where 
both q 2 and p 2 are large ("double tagged"), i.e. 4.5 GeV 2 < p 2 , q 2 < 50 GeV 2 , which only 
gives a very small contribution to the background. Two-photon production of e + e~/ + /~ is 
simulated by the VERMASEREN |TB[ Monte Carlo generator. The different Monte Carlo 
samples are added to provide a complete two-photon sample without double counting. 

Four-fermion processes like conversion and bremsstrahlung diagrams, except for the 
multiperipheral (two-photon) processes, are studied using grc4f. As the studied process 
e + e~ — > e + e~Z/7* is also contained in this class a signal definition is applied (Section ||) 
to classify events either as signal or background. 

Multi-hadronic background e + e~ — > qq is simulated using PYTHIA. As a cross-check 



sample YFS3FF [1/] has been used. Other background processes involving two fermions 



in the final state are evaluated using KORALZ [pl| for e + e — ► /i + /z and e + e — > t + t 



and BHWIDE JTJ] and TEEGG J20| for e+e" -> e+e" (7). 

The integrated luminosity of each of these samples corresponds to at least 5 times the 
data luminosity, except for the two-photon samples, which correspond to at least the same 
as the data luminosity. All Monte Carlo samples are passed through the OPAL detector 
simulation |21j and were subjected to the same reconstruction code as the data. 



The contribution from processes leading to an (e)eqq final state fulfilling the kinematic 
cuts of the signal definition but stemming from other diagrams than the ones shown in 
Figure [I] has been calculated. This (e)eqq background includes processes from multipe- 
ripheral and conversion diagrams shown in Figure |[ For tagged two-photon events, the 
cross-section within our kinematic limits predicted by the HERWIG Monte Carlo simu- 
lation is a = (0.88 ± 0.03) pb. For the conversion processes grc4f predicts a cross-section 
of a — (0.23 ± 0.02) pb within the defined kinematic region. 

4 The momentum transfer squared, q 2 = —Q 2 , in two-photon processes is by definition identical to t 
in our signal process and p 2 = —P 2 is identical to our p 2 . 



4 Event preselection 



The preselection is designed to extract events with two jets and one isolated electron. 
Only tracks and clusters which satisfy standard quality criteria are considered. An algo- 
rithm [^] which corrects for double counting of energy between tracks and calorimeter 



clusters has been used to determine the missing energy and momentum. 

• From the hadronic Z/7* decay two jets are expected in the signal events. For that 
reason the sum of tracks and electromagnetic calorimeter clusters unassociated to 
tracks is required to be greater than 5. 

• All the tracks in the event with an associated electromagnetic cluster of energy more 
than 1 GeV are considered as electron candidates. The ratio of cluster energy to 
track momentum is required to fulfil E e /p c > 0.7. The specific energy loss dE/dx of 
the track in the jet chamber must be consistent with the one for electrons. Rejection 
of electrons originating from photon conversions is implemented using the output 



of a dedicated artificial neural network [E3| . As an isolation criterion no additional 



track in a cone of 0.25 rad half opening angle around the candidate electron track 
is allowed. After subtracting the energy of the candidate the energy deposit in this 
cone must be less than 10 GeV. If more than one electron candidate track satisfies 
these criteria, the one with the smallest additional energy deposit within the cone 
is selected. 

The charge of the candidate has to be consistent with the direction of the missing 
momentum i.e. Q e cos^ m i SS > 0, where Q e is the charge of the track considered as 
an electron candidate and (9 miss the polar angle of the missing momentum. 

Following the signal definition it is required that the invariant mass squared of 
the hadronic system is larger than 25 GeV 2 . The energies and momenta of the 
two jets are obtained from a kinematic fit. The kinematics of the event has to be 
consistent with a topology of two jets and two electrons, with one of the electrons 
going unobserved along the beam pipe. The reconstruction of the jets is performed 
by the k± "Durham" |24|] jet-finding algorithm. The four- vector of the unobserved 



electron is assumed to be (0, 0,p zxmoha , E unohs ), with |p unob s| = ^unobs- Energy and 
momentum conservation are used in the fit within the experimental errors of the 
two jets and the signal electron candidate. An error of 10 mrad has been assigned 
to the direction of the momentum of the untagged electron. The probability of the 
kinematic fit has to be larger than 10~ 6 . 

The e7 centre of mass energy VJ is obtained from the fitted energies and momenta 
of the two jets and the isolated electron. In the signal definition there is an effective 
cutoff on VS at 22.9 GeV. Taking into account the resolution in y§, a cut y§ > 
25 GeV is applied. 

Following the signal definition \t\ has to be greater than 500 GeV 2 . A cut \t\ > 
500 GeV 2 , where t is calculated according to equation ||, is applied. 
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In order to reject Bhabha events the contribution of the three highest energetic 
electromagnetic clusters to the total observed electromagnetic energy is required to 
be less than 93%. Furthermore at least one track with a specific energy loss dE /dx 



in the central tracking chamber |25j not being consistent with an electron hypothesis 
is required. 

Events stemming from interactions of a beam electron with the residual gas or with 
the wall of the beam pipe are not included in the Monte Carlo simulation and are 
rejected by the requirement that the event vertex lies within a cylinder defined by 
I ^vertex I = 10 cm and having a radius of 3 cm. 



After the preselection, 363 events remain in the sample while 339.8 ± 6.9 events are 
predicted by Monte Carlo simulation. The contribution of the signal as predicted by grc4f 
is 84.2 ± 1.7 events. The errors are statistical. 



The overall signal efficiency of the preselection predicted by grc4f is (27.0 ± 0.5)%. 
Splitting up the events into two different kinematic regions defined by the invariant mass 
m qq reveals a dependence of the efficiency on the event topology. In the low mass region 
with an invariant mass m qq between 5 GeV and 60 GeV the efficiency is (21.5 ± 



0. 



Here the main loss in efficiency is due to the multiplicity cut. For m qq > 60 GeV the 
efficiency is (33.8 ±0.7)%. Using the PYTHIA generator similar efficiencies are observed. 

The efficiencies and especially the differences in the efficiencies in the two different 
mass regions can be understood by looking more closely at the topology. The Z/7* is 
predominantly scattered in the forward direction. Therefore in the low invariant mass 
region (7*ee) many particles from the hadronic decay stay in the beam pipe, leading to 
the loss in efficiency due to the multiplicity cut. On the other hand the high invariant mass 
region (Zee) is not affected, as the decay products gain enough transverse momentum to 
be detected within the detector. The differential distribution of the scattering angle of the 
electron is peaked in the backward direction especially for the high invariant mass region, 
strongly reducing the acceptance of the electron. Consequently requiring one electron 
to be detected in the central jet chamber rejects many signal events. The geometrically 
accepted region for the outgoing electron is determined by the minimum number of hits 
required in the jet chamber corresponding to | cos# e | < 0.963. 

The measured distributions of m qq and E e after the preselection are compared to the 
Monte Carlo expectations in Figure |5|. The distributions are well described by the Monte 
Carlo simulation. The resolution of m qq obtained from the kinematic fit is about 3 GeV. 

The main contribution to the background in the low mass region comes from electrons 
from photon conversions in two-photon events. In the high mass region, processes with 
an electron from semi-leptonic W pair decays dominate. In both regions there is a 
contribution from falsely identified electrons. In the signal processes the selected electron 
candidate is almost always the scattered beam electron. Also for tagged two-photon and 
other four-fermion processes mostly correctly identified electrons are found. The tagged 
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two-photon events often satisfy our kinematic signal definition and are difficult to separate 
from the signal process. 

The preselection has been improved with respect to the one applied in Jl[] by including a 
neural network to identify photon conversions, lowering the minimum energy requirement 
for the electron to 1 GeV and a changed isolation criterion. The implementation of the 
kinematic fit improves the resolution in the quantities describing the hadronic system, 
leading to a better description of the kinematic variables m q q, \/§, u and cos#*. 



5 Selection of the signal 



After the preselection, the ratio of signal to background is approximately 1 to 3. In order 
to further reduce the background, the following cuts are applied. The distributions of the 
variables used in each cut are shown in Figure |6|. The numbers of events after each cut 
are shown in Table for data, Monte Carlo signal and the various background processes. 



Cut 


Number of expected events from MC 
Z/7*ee 4f 77 qq 2f Sum 


OPAL 

data 


Presel. 


84.2 ±1.7 


100.0 ±1.8 


106.9 ±6.3 


44.9 ±1.2 


3.8 ±0.6 


339.8 ±6.9 


363 


Cutl 


71.2 ±1.6 


39.9 ±1.2 


80.9 ±4.2 


23.8 ±0.9 


1.2 ±0.3 


217.1 ±4.7 


224 


Cut2 


59.1 ±1.4 


29.5 ±1.0 


54.2 ±3.2 


3.3 ±0.3 


1.2 ±0.3 


147.4 ±3.7 


154 


Cut3 


57.6 ±1.4 


14.2 ±0.7 


53.6 ±3.2 


3.2 ±0.3 


0.7 ±0.2 


129.3 ±3.6 


140 


Cut4 


53.8 ±1.4 


11.4 ±0.6 


35.0 ±2.5 


2.7 ±0.3 


0.4 ±0.1 


103.2 ±2.9 


101 


Cut5 


48.1 ±1.3 


8.3 ±0.5 


9.0 ±1.3 


2.4 ±0.3 


0.3 ±0.1 


68.1 ±1.9 


70 



Table 2: Numbers of expected and observed events for an integrated luminosity of 
174.7 pb~ after each cut. The number of expected signal events is obtained using the 
grc4f generator. The numbers of background events are evaluated using the Monte Carlo 
samples described in the text. The errors are statistical only. 



(Cut 1) The absolute value of the missing momentum must fulfil |p m i SS | > 35 GeV. In 
the signal events the missing momentum is due to the electron which emitted the 
quasi-real photon and remains in the beam pipe. 

(Cut 2) To reduce the background from multi-hadronic events the isolation criterion for 
the signal electron is tightened, requiring that the angle between the electron and 
the second closest track be at least 0.6 rad. 

(Cut 3) For the signal the missing momentum points in the direction of the electron 
staying inside the beam pipe, for this reason the missing momentum vector of the 
event must satisfy QeCos^miss > 0.95. 
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(Cut 4) The maximum energy allowed in the forward detectors Efwd is 30 GeV. With 
this cut the remaining events from the two-photon process where one electron is 
tagged by the forward detectors are reduced. 

(Cut 5) In order to remove the remaining background from tagged two-photon reactions 
those events are rejected where the scattering angle of the electron is in the forward 
direction by requiring — Q e cos# c < 0.65 or E c < 0.35E. This cut is illustrated in 
Figure [7]. 

After all cuts, 70 events are selected while 68.1 ± 1.9 events are expected from the 
Monte Carlo prediction, of which 48.1±1.3 are signal. This corresponds to an overall signal 
efficiency of (15.2 ± 0.4)% according to grc4f. The main contribution to the background 
stems from tagged two-photon events where one of the scattered electrons is detected 
within the detector. In the region removed by the last cut the events from two-photon 
processes are dominant, as shown in Figure |7|. For the region outside this cut, the Z/7*ee 
events are dominant, but still a non negligible contribution from tagged two-photon events 
remains. 

The cross-section a ee is measured in two different regions of m qq , in the low invariant 
mass region 7*ee and in the high invariant mass region Zee. By separating the two mass 
regions at a point where the measured cross-section is near its minimum, the expected 
feed-through, i.e. the number of events with a true value of m qq outside the region it is 
measured in, is very small. The results of the cross-section measurement for both grc4f 
and PYTHIA are summarised in Table |3|. The efficiency for the high mass region is about 
50% larger than that for the low mass region and the expected number of signal events is 
similar. In the low mass region the background is higher, stemming mainly from tagged 
two-photon events. 





grc4f 


PYTHIA 




7*ee 


Zee 


7*ee 


Zee 


Efficiency in % 


13.2 ±0.5 


18.0 ±0.6 


14.4 ±0.7 


19.4 ±0.8 


Expected signal 
Expected background 
Feed through in m qq 
OPAL data 


22.7 ±0.9 

12.1 ±1.1 

0.1 ±0.1 

40 


25.4 ±0.9 

7.9 ±0.8 

0.2 ±0.1 

30 


26.0 ±1.5 

12.1 ±1.1 
0.1 ±0.1 

40 


30.2 ±1.6 

7.9 ±0.8 

0.3 ±0.2 

30 


Measured cross-section in pb 


1.20 ±0.28 


0.69 ±0.18 


1.11 ±0.26 


0.64 ±0.17 


Predicted cross-section in pb 


0.98 ±0.01 


0.80 ±0.01 


1.03 ±0.02 


0.88 ±0.02 



Table 3: Comparison of OPAL data with the cross-section a ee predicted by the grc4f and 
the PYTHIA Monte Carlo generator. The errors are statistical only. 



In Figures |] and |9] the measured event distributions for several variables are com- 
pared to the ones predicted from the grc4f Monte Carlo simulations. For comparison the 



event distributions from PYTHIA are given in Figure |10| for some variables. The limited 
statistics of the data does not allow to distinguish between the two simulations. 
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6 Determination of differential cross-sections 



To compare the results of this analysis with results from other experiments, differential 
cross-sections are calculated from the distributions of the event variables. For the observed 
process differential cross sections dcr ee and dcr e7 have been determined. 



6.1 Differential cross-sections der, 



ee 



For the determination of the differential cross-sections the experimental resolution is of 
importance. If the experimental resolution is much smaller than the chosen bin width 
and the distribution is flat, then the correlation matrix will be close to the unit matrix. 
But if the distribution is peaked, like for m qq - then a large fraction of events measured in 
bins around the peak originated from bins other than the one they had been generated 
in. Therefore the correlation for each variable between the generated and the measured 
distribution has to be determined. 

The correlation matrix M between the generated and the measured distribution has 
been calculated for each variable shown in Figures |8| and ^ using the bin width shown 
there. The matrix M is given by: 

"«•'"> =iH6j (n) 

and fulfils the following condition: 

N gen (t) = Y / M(i,j)N det (j), (12) 

3 

where N gen (i) is the number of events generated in bin i, Nd e t{j) is the number of events 
measured in bin j and G(i,j) is the number of events generated in bin i and measured 
in bin j. The matrix M has been determined using the grc4f MC. For most variables M 
is very similar to the unit matrix with some small off-diagonal elements. For m qq around 
the Z-mass non zero elements exist also away from the first off-diagonal. 

The matrix M has also been calculated from the PYTHIA signal MC and no significant 
difference with the one determined from grc4f has been observed. 

The differential cross sections are then given by 

^ = £M(z,j)(AU.7>) -^(*))^Z^l(i) ' (13) 

where x is the variable used, C ee the integrated luminosity and e(x) the efficiency in a 
given x bin. The differential cross-sections are shown in Figure [U]. 
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6.2 Differential cross-sections d<x, 



e7 



A further aim of this analysis is to calculate the differential cross-sections d<r e7 . This 
allows the results from a given e + e~ centre-of-mass energy to be compared with other 
energies as well as with results from other colliders. To calculate the differential cross- 
section for a variable x, Equation [I] has to be inverted: 



da, 



c 7 



dx 



1 


1 


da ee 


Az 


D ei [z, u) EPA 


dx 



(14) 



with D ei (z,u) E -p- A being the photon-flux and Az = As/s. The lower limit of y§ is given 
by the signal definition as 23 GeV and the upper limit is 160 GeV, resulting in Az 
To calculate d<7 e7 /dx the mean of the inverse of the photon-flux < 1/D, 



z,u 



EPA 



0.702. 
> is 



[z,u 



> 



calculated in bins of u and s. Taking into account the dependence of < l/D^y*, u JEPA 
on u is necessary, as the efficiency varies with u. We have chosen three bins in u and 
four in s. The bin boundaries as well as the mean values of < 1/D e ^(z, u)^^ > ar e given 
in Table f|. The dependence of < 1/D e ^(z, u) E -px > on u is small while it is large for s. 
The width of the bins is chosen to be at least three times larger than the experimental 
resolution. 



/ 1 \ 


bins in u in GeV 2 




\D e ^{z, Mj E p A / 


-50000 to -1000 


-1000 to -200 


-200 to -10 


bins 23 to 50 


0.69 


0.67 


0.59 


in Vs 50 to 80 


2.02 


2.03 


2.24 


in GeV 80 to 110 


6.08 


6.61 


7.19 


110 to 160 


13.90 


14.49 


16.06 



Table 4: Mean values of the inverse of the photon flux < 1/D e ^(z, w)kpa > ^ n bins of u 
and 



The differential cross-section d<j ee /dx is calculated according to Equation [13| using bins 
of u and s. 



da cc 
dx 



E 

u,s 



N(s,u,x)/Ax 

£ee e(u,x) 



(15) 



where N(s,u,x) is the number of events in bins of s, u and x after subtracting the 
background and using the matrix M. Larger bins compared to the previous section have 
been used and M is calculated using these bin sizes. The efficiency e(u, x) is calculated 
only in bins of u and x as it is flat in Vs. This results in the differential cross-section in 
the e7 system to be given by: 



da, 



c 7 



dx 



E 



u,s 



N(s,u,x)/Ax 1 



Co 



e[u,x\ 



Az 



U e ^{z, u) E p A 



(16) 
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The measured differential cross-sections are shown in Figure O and are compared to the 
generated distributions. 

For calculating the total cross-section a n (\/I) as a function of yl the same method 
as above is applied, with the difference that Az is calculated for each bin of Vs. 



7 Systematic error studies 



For the calculation of the total and differential cross-sections the efficiencies, unfolding 
matrix and backgrounds are taken from Monte Carlo simulations. It is therefore important 
to study systematic effects resulting from these simulations. 



7.1 Systematic error studies for the total cross-section 

The systematic errors on the efficiencies come mainly from imperfect modelling of the 
detector response. This can lead to discrepancies between the data and the Monte Carlo 
simulation in the distributions of the cut variables. These systematic errors can be esti- 
mated, for example, by comparing Monte Carlo simulation and data. 

This has been done using the events selected by the preselection. For these events each 
of the selection cuts has been applied separately and the relative difference in the number 
of events selected in data and in Monte Carlo has been assigned as a systematic error 
after quadratically subtracting the statistical error. In cases where this results in a value 
being lower than the statistical error, conservatively the statistical error is used. An error 
common for the whole mass range has been calculated and no distinction between the low 
and high mass range has been made. The values of the errors are listed in Table [|. These 
systematic uncertainties are used for both the 7*ee-like and Zee-like kinematic regions 
since the efficiencies of each of these selection cuts are similar in the two regions. 

As a cross check the systematic errors have also been estimated by comparing Monte 
Carlo simulation and data for the process W + W~ — ► qqe^ e which has the same observable 
final state as the Z/7*ee process. W + W _ — > qqez/ e events are selected according to the 
procedure described in [^6| and then each selection cut is applied separately to this sample. 
For those cuts where the distribution of the cut variables is similar for W + W~ — ► qqez/ e , 
Zee and 7*ee events (the absolute value of the missing momentum, the electron isola- 
tion and the electron angle and energy) no difference within the statistical error in the 
systematic error compared to the method described above has been observed. 

The relatively smaller overall efficiency for 7*ee-like events arises mainly from the 
multiplicity cut in the preselection. To asses the systematic uncertainty of this cut the 
number of tracks and clusters required has been changed by ±1. 
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7*ee 


Zee 


multiplicity 


0.024 


— 


bmissl > 35 GeV 


0.042 


0.042 


electron isolation 


0.038 


0.038 


— QeCOS^e CUt 


0.028 


0.028 


QeCOS^miss > 0.95 


0.033 


0.033 


^fwd < 30 GeV 


0.056 


0.056 


detector simulation 


0.094 


0.091 


efficiency 


0.049 


0.050 


background 


0.038 


0.033 


Total 


0.113 


0.109 



Table 5: Relative systematic uncertainties of the cross-section measurements. The entry 
"detector simulation" is the quadratic sum of the signal efficiency uncertainties for the 
single cuts listed in the rows above it. 

The uncertainty in the efficiency due to the choice of a particular Monte Carlo event 
generator is estimated by comparing PYTHIA and grc4f. For this comparison only events 
with the signal electron within the acceptance of the detector, defined by a cut on gener- 
ator level on the angle of the electron | cos6^ cn | < 0.963 are used. The relative difference 
in efficiency of the two different Monte Carlo generators after subtracting the statistical 
errors quadratically is taken as a systematic error. 

Uncertainties affecting the residual background have been evaluated separately for 
each of the three main background classes remaining after all cuts. Background-enriched 
samples are obtained by inverting or omitting one or two cuts, while the other cuts remain 
unchanged. The full difference between the number of events remaining in the data and 
the number of expected events from Monte Carlo is taken as a systematic uncertainty. 
For the background from four-fermion final states the cut on the fit probability is omitted 
and the cut on the electron isolation is inverted. After applying these cuts, a relative 
difference of 17% between the data and Monte Carlo is observed. The background from 
tagged two-photon events is enriched by inverting the cut on the electron's angle and 
energy. The relative systematic error is 10%. By omitting the cut on the angle of the 
missing momentum and inverting the electron isolation cut the multi-hadronic background 
is enriched, leading to a relative systematic uncertainty of 12%. 

The numbers of background events after all cuts in the 7*ee region are 4.47±0.39±0.75 
from four fermion events, 1.07 ± 0.19 ± 0.11 from multi-hadronic processes and 5.80 ± 
1.01 ±0.73 events from tagged two-photon reactions. In the Zee region the corresponding 
contributions are 3.81 ± 0.34 ± 0.64, 1.38 ± 0.22 ± 0.14 and 2.46 ± 0.66 ± 0.31 events. 
Additional background sources contribute with less than 1 event for each mass region. 
This leads to a relative error on the cross- section of 3.8% in the low mass region and 3.3% 
in the high mass region, as quoted in Table [5|. 
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In further studies contributions to the background not modelled in the Monte Carlo 
were investigated. Using random beam-crossing events, the interactions between the beam 
particles and the gas inside the vacuum pipe were found to be negligible. 

For the multi-hadronic background an additional systematic cross-check is applied by 
comparing the prediction of two different Monte Carlo generators. A good agreement 
between the PYTHIA and YFS3FF Monte Carlo generators has been found within the 
statistical errors. Consequently no additional systematic error has been assigned. 



7.2 Systematic error studies for the differential cross-sections 

The systematic errors for the differential cross-sections stem mainly from the imperfect 
detector simulation in the Monte Carlo and from the uncertainty in the unfolding matrix 

M. 

The error assigned for the detector simulation is taken to be the same as for the total 
cross-section. A value of 9.4% is assigned. The unfolding matrix M is calculated using 
both the grc4f and the PYTHIA Monte Carlo sample and the relative difference between 
the two is taken as a systematic error. For the background the same errors as determined 
for the total cross-section are used. The errors for the three different background classes 
are taken into account in each bin of the event distributions. The efficiency e is calculated 
in bins of the variable x for the differential cross-sections da ee and in bins of u and the 
variable x for the differential cross-sections d<T e7 . It has therefore much larger statistical 
errors than in the calculation of the total cross section. Within the statistical errors no 
difference between the efficiencies determined from the grc4f and the PYTHIA Monte 
Carlo samples has been observed. Consequently no systematic error has been assigned. 

Due to the small statistics of the data sample the systematic error for the differential 
cross-sections is much smaller than the statistical one. 



8 Results and discussion 



The cross-section for the process e + e~ — > e + e~Z/7* has been measured in a restricted 
phase space, defined by Lorentz invariant variables, in two different regions of the mass of 
the hadronic system, corresponding to either a 7* or a Z° in the final state. With the cut 
at m q q = 60 GeV the 7* and the Z° are well separated. With an integrated luminosity 
of about 175 pb -1 , a total of 70 candidate events have been observed, while 20.0 ± 1.4 
background events and 48.1 ± 1.3 signal events are predicted, giving a total of 68.1 ± 1.9 
events. The cross-sections times branching ratio for the decay of Z/7* into hadrons, <T ee , 
are found to be a = (1.20 ± 0.28 ± 0.14) pb for 7 *ee and a = (0.69 ± 0.18 ± 0.08) pb for 
Zee final states within the kinematical definition listed in Table |l]. A large part of this 
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cross-section is not detectable within the detector as the scattering angle of the electron 
is in the very backward direction. For the calculation of the cross-sections the efficiencies 
predicted by the grc4f generator are used. The cross-sections measured using efficiencies 
predicted by PYTHIA lie well within the errors. 

The distributions of E q q, m q qand E e are shown in Figures ||(a) to (c), respectively. 
The distribution of £" q q shows two peaks, one at the beam energy and one at about 
115 GeV. From equation [| one obtains for the largest part of the cross-section at u = 

and QeCOS^qq = — 1 

771 - 

£ q q = -Ebeam + . p ~ • (17) 

Consequently, the peak at the beam energy corresponds to 7*ee and the peak at 115 GeV 
to the Zee process. The tail at lower energies is due to the parts of the cross-section 
where u^O and Q e cos# q q ^ — 1. Two peaks are visible in the invariant mass distribution, 
stemming from the contributions of the two gauge bosons, the 7* and the Z. The expected 
background is flat over the whole mass range up to 100 GeV. For both distributions the 
shapes of Monte Carlo and data are in good agreement. The two Monte Carlo generators 
give similar distributions and more statistics is needed to distinguish between the two. 

In Figure |9| the distributions of the scattering angle 6* of the Z/7* in the e7 rest-system 
with respect to the incoming photon direction, as well as the distribution of the kinematic 
invariants \/§, i and u are shown and are compared with the predictions of grc4f. The 
scattering angle 9* peaks strongly in the backward direction and the agreement between 
data and Monte Carlo is good. The observed structure of the distribution of v«, can be 
understood in terms of the 7*ee final states in the low V§ region and the Zee in the high 
Vs region. In Figure [10] the measured distributions are compared to the distributions 
from PYTHIA. The predictions of both grc4f and PYTHIA are in agreement with the 
data. The distribution of t is peaked towards and shows a long tail towards large values. 
The distribution of u shows the typical behaviour of a w-channel process, a peak at zero. 

From these event distributions the differential cross-sections dcr ee are derived and are 
shown in Figure [11]. Only Lorentz invariant quantities have been derived. The Monte 
Carlo describes the data well, except for small values of Vs, where the Monte Carlo 
underestimates the data. In the distribution of dcr ee /dm q q for small values of m q q the 
steep falloff of the cross-section with increasing invariant mass m q q as well the peak at 
the Z-mass are very well visible. The distribution of d<7 ee /dVs shows a decrease with 
increasing Vs until the threshold for Z Boson production is reached. 

From the differential cross-sections dcr ee the differential cross-sections for the sub- 
process e7 — > e Z/7* are deduced, based on a factorisation ansatz using the modified 
Equivalent Photon approximation from Equation [5] They are shown in Figure |T2 . 



The distribution of d<Je 7 /dm qq shows a strong peak around the Z Boson mass. The 
differential cross-section da ei /dit shows a very sharp peak at 0, as expected for this u- 
channel process. For d<r e7 /d£ the Monte Carlo does not describe the data so well. The 
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increase for t towards is well reproduced, but for large negative values of t the Monte 
Carlo lies constantly above the data. 

In Figure 0(b) the total cross-section a ei {\fs) is shown. Note the dip in the cross- 
section at \/I around 70 GeV. The total cross-section a n (V§) shows a decrease with 
increasing Vs until the threshold for Z Boson production is reached. This is the first 
measurement of <r e7 (Vs) around the Z Boson threshold. a e7 (Vs) is independent of the 
e + e~ centre-of-mass energy and can therefore be compared with measurements at other 
e + e~ centre of mass energies as well with measurements at other colliders, e.g. HERA. 

Within the statistical error, the Monte Carlo predictions are in good agreement with 
the data. But there is a tendency that the data are higher than the Monte Carlo in the 
low v§ region, while they are too low in the high Vs region. 



9 Conclusions 



The process e + e~ — > e + e~Z/7* and its subprocess ej — > eZ/7* have been studied. For 
the process e + e~ — > e + e~Z/ / y* the cross-section times branching ratio for the decay of the 
Z/7* into hadrons at y/s = 189 GeV has been measured within a restricted phase space 
to be a = (1.20 ± 0.28 ± 0.14) pb for m^< 60 GeV and a = (0.69 ± 0.18 ± 0.08) pb 
for m q q> 60 GeV. The Monte Carlo generators grc4f and PYTHIA both predict cross- 
sections within one standard deviation of the measured values. The cross-section a^(vs) 
for the subprocess e7 — > eZ/7* has been determined in a range of V^ from 23 to 160 GeV. 

Differential cross-sections da ee and d<r e7 have been determined and compared to the 
ones from the Monte Carlo generators grc4f and PYTHIA. The generators describe all 
distributions well. 
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Figure 1: Diagrams for the process e + e — > (e)eZ/Y 





Figure 2: Further diagrams leading to the final state eeqq. On the left is the multiperipheral 
diagram and on the right is the conversion diagram. 
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Figure 3: Distribution of the charge weighted cosine of the scattering angle of the electron 
— Q c cos# c versus its energy E c as predicted by the grc4f signal Monte Carlo at generator 
level. The solid line corresponds to \t\ = 500 GeV 2 . 
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Figure 4: Distribution of the transverse momentum p t of the Z/7* in the cy rest-frame, 
the hadronic mass, m q q ; the electron energy, E e , and t at generator level for the signal 
for all events fulfilling the signal definition. The histograms show the distributions for the 
grc4f sample and the points for the PYTHIA sample. Only the statistical errors of the 
PYTHIA sample are shown. The distributions are normalised to the data luminosity. 
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Figure 5: Distribution of (a) m q q and (b) the electron energy E c after the preselection. 
The histograms show the contributions from the various processes and the points represent 
the data. The signal contribution is taken from the grc^f Monte Carlo sample. Only 
statistical errors are shown. 
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Figure 6: Distributions of variables used in cuts in the preselection and in the selection, 
(a) P , the fit probability, (b) p ni i SS , the missing momentum, (c) the angle between the 
electron track and the nearest track, and (d) E FWD , the energy deposited in the forward 
calorimeter before applying the cut on that variable. The arrows indicate the selected 
region. Only statistical errors are shown. The histograms show the contributions from the 
various processes and the points represent the data. 
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Figure 7: Distribution of the cosine of the scattering angle of the electron cos6 c versus its 
energy E e before the last cut. In (a) the distribution is shown for the signal, in (b) for 
the data and in (c) for the tagged two-photon Monte Carlo sample. Both Monte Carlo 
samples have been scaled to the same integrated luminosity as the data. The solid line 
indicates the cut in the preselection at \t\ > 500 GeV . The area inside the dashed line 
corresponds to the region discarded by the last selection cut. 
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-E q q and (b) mass m qq - of the quark system and of (c) the electron energy E e . The open 
histogram shows the signal simulated with the grc4f Monte Carlo, the hatched histogram 
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Only statistical errors are shown. 
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Figure 9: Distribution of (a) the scattering angle in the cy rest-frame cos 6*, and the 
kinematic invariants (b) y§, (c) i and (d) u after all cuts. The open histogram shows the 
signal simulated with the grc4f Monte Carlo, the hatched histogram shows the backgrounds 
from various Monte Carlo simulations and the points the data. Only statistical errors are 
shown. 
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Figure 10: Distribution of (a) the mass m q q and the kinematic invariants (b) vs, (c) i 
and (d) u after all cuts. The open histogram shows the signal simulated with the PYTHIA 
Monte Carlo, the hatched histogram shows the backgrounds from various Monte Carlo 
simulations and the points the data. Only statistical errors are shown. 
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Figure 11: The differential cross- sections (a) da ee /dm q q, (b) da ee /dv§, (c) da ee /di and 
(d) da ee /du are shown. The open histograms show the signal simulated with the grcJ^f 
Monte Carlo generator and the points the data. The errors show the statistical and sys- 
tematic errors added in quadrature. The contribution of the statistical errors is indicated 
by the horizontal bars. 
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Figure 12: The differential cross-sections (a) c& C7 /dm q q, (c) da ei /dt and (d) da cl /du are 
shown. In addition (b) the cross section <t C7 (Vs) is shown The open histograms show the 
signal simulated with the grc4f Monte Carlo generator and the points the data. The errors 
show the statistical and systematic errors added in quadrature. The contribution of the 
statistical errors is indicated by the horizontal bars. 
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